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Waveguide design has always been an important part 
of photonic integrated devices, and its research is also 
developing. In this paper, the coupling between waveg- 
uides is studied from two aspects of coupled mode the- 
ory and finite element simulation. In addition, this pa- 
per gives an example of dispersion engineering caused 
by dual-core waveguides structure, which has a wide 
range of applications in supercontinuum. To sum up, 
coupled mode theory has important reference signifi- 
cance for the design of modern photonic devices. 
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1. INTRODUCTION 


Evnascently light coupling among optical waveguides repre- 
sents one of the most basic and applied effects in photonic cir- 
cuits as well as in fiber optics, and have enabled a diversity 
of fundamental building blocks including directional couplers 
[1, 2], mode converters [3, 4], multiplexers [5, 6], interferometers 
[7], as well as resonators [8, 9], and have allowed for the ob- 
servation of electromagnetically induced transparency [10, 11], 
integrated optical isolation [12], Parity-time symmetry [13], etc. 
In the phase matching condition, the light coupling would give 
rise to a number of distinct mode states as the superposition of 
individual modes, e.g. the pair of symmetric and anti-symmetric 
modes in a classic two-mode system. Such states are analogous 
to discrete energy states in molecular systems, and similar-to- 
molecule performances have been observed and controlled in 
photonics ways [14, 15]. On aspect of nonlinear photonics, such 
states have enabled advanced control of light in terms of the 
dispersion engineering [16, 17] and the nonlinear interaction. 
Photonic supercontinuum generation is another scheme that 
could benefit from the light coupling effect, where localized 
dispersion engineering is enabled by induced dispersion of the 
mode coupling [18-20]. In particular the induced anomalous 


dispersion from the anti-symmetric mode is usually desired for 
soliton-regime supercontinuum, which could be designed at 
arbitrary wavelengths including the mid-infrared. Nevertheless 
the supercontinuum could also model a broadband and unde- 
pleted pumping system, where weak and localized mode cou- 
plings (i.e. coupling with higher-order modes or to higher-order 
harmonics) would feature resonant absorption and leave narrow 
line-shapes in the spectrum [21, 22], which again is analogous to 
molecular absorption usually featured in spectroscopy. 


Recently, a mid-infrared dispersion engineered and flat su- 
percontinuum has been successfully implemented in a dual-core 
coupled waveguide structure, and was further applied in a dual- 
comb spectrometer for parallel gas-phase detection [23]. Yet, the 
design process for such coupled waveguide structure is non- 
trivial. The most applied criterion is the coupled mode theory 
[24], which treats a coupled waveguide element (core) as di- 
electric disturbance, and predicts the light coupling dynamic 
on basis of individual modes (in individual waveguides). In- 
deed, the theory has been well studied and widely applied in the 
design of fundamental optical components, as well as recently 
in topological photonics [25, 26], photonic crystals [27, 28] and 
quantum cascade lasers [29]. Yet, the quantitative accuracy of 
the coupled mode theory remains questionable, especially to- 
wards the mid-infrared region where the mode confinement of 
the waveguide is reduced. As such, the modern way of perform- 
ing the dispersion engineering is usually based on numerical 
simulations by means of e.g. the finite element method, which 
would certainly be more time-costing compared with the analyt- 
ical method as repeated simulations are required to explore the 
waveguide dispersion on each structural parameter. This posts 
limitations to systematic. 


In this letter, we rephrase the coupled mode theory, in which 
the coupling induced self-correction on the propagation con- 
stant is highlighted, and is understood as the key to extend 
the accuracy of the theory towards the cut-off (i.e. to the long 
wavelength range). With nonlinear effects, a revised coupled 
nonlinear Schrödinger equation is also derived. We then show- 
case the engineering of dispersion from the mode coupling, in 
the mid-infrared range in a photonic integrated dual-core waveg- 
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uide, which has excellent accuracy and could be applied in the 
prediction of mid-infrared supercontinuum generation includ- 
ing the dispersive wave generation. 


2. THEORETICAL MODEL 


Here, we limit the coupled mode analysis regarding the most 
typical waveguide structure, namely the structure with two 
parallel waveguide cores. The coupling of mode amplitudes 
(a(z) and b(z)) in individual cores , along the propagation axis 
(z), are generally described as: 

date) = = —ikapbe P 
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where Kap and Kpa are the ame coupling coefficients between 
the modes, and Af indicates the phase mismatch. The complete 
form of Eq. 1 could be derived based on the dielectric perturba- 
tion theory. In this context, the electric field of the general wave 
in the dual-core waveguide can be approximated by: 


E(x,y,zZ,t) = 
alz)Ea(x, yet Eetu] + B(z)Ep(x, yer] (ay 


where E, p(x, y) indicates the transverse mode field distribution. 
It should be noted that a self-correction factor (i.e. Kaa OF Kpp) 
is introduced, which physically is resulted by dielectric pertur- 
bations to one of the modes due to the presence of the other 
one. 

In principle, the electric field given by Eq. 2 must obey the 
wave equation 


3 Pa 3 2 2 
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where the refractive index (ri) distribution is further expanded 
as: n(x, y) = n(x, y)+ Anž(x, y)+ Ani(x, y) (ns is defined as the 
ri of me substrate or cladding area outside the waveguide cores, 
and An? ab indicates the difference of squared ri of the core to that 
of the substrate. Thus, Eq. 3 can be expanded, and a dynamic 
with respect to the mode amplitude a(z) is obtained as: (similar 
for the dynamic of b(z)) 


a = —ia(Kaa a 2 ff An;( (x, y)E,E, dxdy) 
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where €9 is vacuum permittivity. 

To force a consistent form with Eq. 1, one could conclude the 
coupling coefficients, the self-correction factors, and the phase 
mismatch as the following: 


Kab = on Jf And(x,y)EpE* dxdy 

Kba = a ff An +(x, y)Ea E,dxdy 

Kaa = 2 ff An a(x, y)Ea E,dxdy (5) 
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We further make A = Kpaae "P aKa and B = Jragbe™? oR) 
So a set of coupled mode equations is derived as: 
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where the renewed coupling coefficient K = YKabKba naturally 
form a complex conjugate pair in the transfer matrix. As such 
the energy conservation is ensured for the coupling between A 
and B. 

Moreover, the eigenvalue(s) of the transfer matrix could be 
calculated, i.e.: 


(Ba — Kaa) + (Bo — Kob) 


À12 = 7 


S 
+ z (7) 


where S = yA? + 4x?. The eigenvectors are: 
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Physically, such eigenvalues are comparative with the propa- 
gation constants of the symmetric and anti-symmetric modes in 
the dual-core waveguide, and the elements in eigenvectors indi- 
cates the amplitude of individual modes (A and B) in forming 
such superpositioned modes (supermodes). 

To facilitate the simulation of nonlinear pulse propagation 
and supercontinuum generation in such a coupled dual-core 
waveguide system, Eq. 6 is further extended to contain the 
nonlinear phase modulation effects, forming a set of coupled 
nonlinear Schrödinger equations (NLSEs), i.e.: 


HE = —i(Ba — Kaa — Bs)A’ - ixB' - iyaF [|A'Y A'] (9) 
dE = —i(By — Kob — Bs)B' - ix A’ — ino [|B'} B'] 


where a dispersionless phase constant f; is assumed underlying 
the mode amplitude, i.e. A' = Ae and B' = Be fe, Bs = 
Ba(ws) + (w - Ws)0g" (ws is the central angular frequency of the 
pumping pulse, and vg is a referenced group velocity of the 
pulse), Yap indicates the nonlinear coefficient for the self-phase 
modulation, which is dependent on the effective mode area 
(Aeff a,b) in individual waveguide cores, i.e. Yq, & Ace ae 

Compared with the conventional NLSE, a primary change 
of Eq. 9 is the presence of the self-correction coefficients, which 
may impact the dispersion profile and lead to deviations on 
the nonlinear pulse propagation dynamic, including the pulse 
seeding supercontinuum process. 


3. RESULTS 


A. Symmetrical double waveguides 


For dual-core waveguide structure as shown in Fig1.(a) (h=1000 
nm, w=1000 nm, g=1000 nm; nab = Nsi,N,, Ms = Nsio,). When 
the same waveguides are closer, the optical mode propagating 
in one will be coupled to the other, which effectively changes its 
phase i.e. the propagation constant(A(A), A is the wavelength of 
the light). Physically, mode coupling will bring the hybridization 
of mode-field distributions, resulting in a pair of supermodes, 
namely the symmetric and anti-symmetric modes. 

Usually, the propagation constants of super modes fs1,s2 can 
be calculated by COMSOL simulation. However, the eigenval- 
ues A; solved by the Eq. 7 can also be used as the propagation 
constant of super modes. In order to study the influence of cou- 
pling, uniformly observe the change of propagation constants. 
Through data processing, it is found that the contrast satisfies 
the following relationship Bs1 — By = A1 — Bp, Bs2 — Ba = A2 - Bay 
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Fig. 1. Coupling between the same waveguides. (a) Design of 
dual-core waveguide structure: h is height, w is width, and g 
is the gap between two waveguides. (b) Comparison of propa- 
gation constant changes; Inset is the electric filed distribution 
of super modes s1 and s2. (c) Error caused by changing the 
gap. The blue line and red line represent 2500 nm and 3000 nm 
wavelength positions respectfully. 


which can be clearly observed in Fig1.(b). In addition, the ob- 
tained self coupling coefficient is very small whose effect on 
the propagation constant usually be neglected. In fact, if do not 
make the ba — Kaa for correction as mentioned in Eq. 2, A1 — Bp 
will directly equal to Kap which will bring the inaccuracy as 
shown in Fig1.(b). 

By changing the gap between two waveguides to smaller, 
as shown in Fig1.(c), it is found that the relative difference of 
the relationship will gradually increase. Especially, the asym- 
metry results of two supermodes in this structure will become 
more and more obvious at the long wavelength, which can be 
explained by the existing of self coupling coefficient. 


B. Asymmetric double waveguides 

Since different waveguides with different propagation constants, 
2 suitable waveguides that exist the coupling should be selected. 
As shown in Fig2.(a), the start wavelength of the coupling and 
the phase matching point of the dual-core waveguide structure 
can be easily found by calculating the relative effective refractive 
index(An; = nj; — na). Finally the designed waveguide a and 
b can be determined(width of a is 1000 nm, b is 3000 nm and 
g=1000 nm, h=1000 nm). 

At the same time, described as A\ng2 in Fig2.(a), the appear- 
ance of anomalous dispersion caused by coupling can be easily 
observed. In addition, Fig2.(b) shows the variation value of 2 
independent electric field composition which is related to the 
eigenvector mentioned in Eq. 8. For example E; = m-E, +n- Ep, 
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Fig. 2. Coupling between 2 different waveguides. (a) Rela- 
tive effective refractive index difference of super modes s1,s2 
and single waveguide a,b calculated based on waveguide a; 
(b) Coefficients obtained from eigenvectors for synthesizing 
supermodel electric fields. (c) Dispersion calculation of inde- 
pendent waveguide a and dual-core waveguide structure with 
different gaps. 


where |m| + |n| = 1. Especially, when two waveguides are the 
same, |m| = |n|. Besides, for dispersion engineering, changing 
the gap can adjust the region of anomalous dispersion as shown 
in Fig2.(c). 

For supercontinuum generation, Fig3.(a) shows the change of 
relative propagation constant of 2 independent waveguides and 
dual-core waveguide structure which are calculated by COM- 
SOL simulation. In Fig3.(b), pumping waveguide a at 2000 nm, 
the supercontinuum spectrum after coupling can be achieved 
from Eq. 9. Besides, the comparison between simulation and 
theoretical calculation of supercontinuum generation also be 
given in Fig3.(c). Especially, the grey line gives the result of 
without self-correction factor (i.e. Kaa OF Kpp). 

According to the supercontinuum, although there are some 
differences, the correction of introducing the self-correction fac- 
tor is necessary, which can be observed especially from the posi- 
tion of the dispersive wave. Importantly, we can draw a conclu- 
sion that when designing waveguides based on coupled mode 
theory will not only improve the accuracy but also speed up the 
calculation. 


4. CONCLUSION 


From the comparison results, we can get some meaningful con- 
clusions. First, since waveguide modes tend to be cut off at 
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Fig. 3. Supercontinuum generation. (a) Calculated dispersion 
landscape of the anti-symmetric mode of dual-core waveguide 
structure (gap = 600 nm) and the fundamental mode of waveg- 
uide a,b. (b) Supercontinuum generated by two independent 
waveguides based on coupled mode theory. (c) Comparison 
of supercontinuum generation based on coupled mode theory 
and software simulation. 


long wavelength, using the coupled mode theory can effectively 
avoid the errors caused by simulation; Second, the inconsistency 
of the calculated supermodes for the dual-core waveguide struc- 
ture of the same waveguides can be explained by the existence 
of self coupling coefficient. 

Last but not least, we focus on the coupling between 2 dif- 
ferent waveguides which can bring abnormal dispersion for 
dispersion engineering. And we find the differences between 
theoretical calculation and simulation calculation, especially the 
effect of self coupling coefficient. Nevertheless, we know the 
difference between the supercontinuum spectrum generated by 
the two methods, and we can use it for reference in the future 
waveguide design work. 
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